Superconductivity-induced macroscopic resonant tunneling 
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We show analyticaUy and by numerical simulations that the conductance through 7r-biased chaotic 
Josephson junctions is enhanced by several orders of magnitude in the short-wavelength regime. We 
identify the mechanism behind this effect as macroscopic resonant tunneling through a macroscopic 
number of low-energy quasi-degenerate Andreev levels. 
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Resonant tunneling is the process by which quantum 
tunneling is enhanced by resonant transfer through inter- 
mediate quasi-bound states [l[ . The paradigmatic exam- 
ple is a one-dimensional double-barrier structure, where 
the transmission probability is given by 

T1T2 



I + R1R2- 2y/lhThcos[Q{E)] ' 



(1) 



in terms of the transmission and reflection probabilities 
Ti 2 = 1 — i?i,2 of the individual barriers, fn the tun- 
neling regime, T12 <C 1, narrow quasi-bound states exist 
between the two barriers, with well resolved quantized 
energies, em- When the energy of the tunneling particle 
coincides with one of these energies, Q{E = e,„) = 0, and 
in the case of symmetric barriers, Ti ~ T2, the transmis- 
sion is perfect, Tj^csi^m) = 1- This is to be contrasted 
with the transmission probability T{E) = T1T2/4 away 
from resonance, and the incoherent transmission proba- 
biUty T{E) = TiT2/(Ti +T2) one obtains when inelastic 
scattering occurs between the two barriers Q • 

Resonant tunneling also occurs in higher dimensions. 
In chaotic systems with no spatial symmetry, there is no 
degeneracy of the intermediate states. Therefore, con- 
sidering linear transport at low temperature, resonance 
occurs with at most one intermediate state at a time, 
leading at best to an increase of the average conductance 
by an amount Go = 2e^ jh - it is a microscopic effect of 
order one. In this article, we show that the proximity 
of the intermediate system to two superconductors can 
lead to a totally different phenomenology, where resonant 
tunneling through a macroscopic number oc of inter- 
mediate levels occurs at the Fermi energy. This results 
in a conductance G oc GqN^ at resonance, much larger 
than the nonresonant conductance cx GoTnNn- The reso- 
nance condition is met when the phase difference between 
the two superconductors is = tt. We foresee that this 
macroscopic resonant tunneling effect might have appli- 
cations in current switching devices and magnetic flux 
"transistors" . 

The system we investigate is sketched in Fig. [T] A bal- 
listic metallic quantum dot is connected to two metallic 
electrodes (L and R, each carrying 1 channels) and 

two superconducting electrodes {Si and S2, each carrying 
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Figure 1: (Color online) a) Andreev reflection by a su- 
perconductor (S). An incoming electron (e) in a normal 
metal (N) near the Fermi energy Ef is reflected as a hole 
(h) with opposite velocity, b) Schematic of our system. A 
ballistic metallic cavity is connected to two normal (L and 
R) metallic leads, and two superconducting {Si and ^2) 
leads. The superconducting pair potentials have a phase 
difference of (p. Two periodic Andreev orbits contributing to 
macroscopic resonant tunneling at = tt are drawn. 



A^s ^ 1 channels). The electrodes are coupled to the dot 
via tunnel contacts of transparency < T^^g ^ Ij such 
that 1 < FnA^n < A^n, and FgTVs > 1. We are interested 
in transport between the two normal leads as a function 
of the phase difference (f> between the two superconduc- 
tor's pair potentials, A^^ = cxp[—i(j>], \Asi 2I = A. 

The physics in our system is to a large extent governed 
by Andreev reflection [3|] . At low energy, this is the dom- 
inant reflection process at an interface between a metal 
and a superconductor, where an electron is retroreflected 
into a hole, and vice-versa. The process is sketched in 
Fig. [TJ When the excitation energy e is negligible against 
A and the Fermi energy Ep, the retroreflection is perfect 
and the hole (electron) exactly retraces the path previ- 
ously followed by the electron (hole), with an additional 
Andreev reflection phase of — 7r/2 ± 0/2. In absence of 
normal lead (one then has an Andreev billiard [J), An- 
dreev reflection renders all classical paths periodic in a 
cavity that would be chaotic otherwise. 

When the cavity is weakly connected to external leads, 
transport can be resonantly mediated through those peri- 
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odic orbits that touch the contacts to the leads and both 
superconductors. Two such orbits are depicted in Fig. [1] 
Each such orbit represents a family of scattering trajecto- 
ries constructed from a primitive trajectory, and an An- 
dreev loop that can be traveled p times, p = 0,1,2,.... 
We indicate segments of trajectories as 7a"^ with a su- 
perscript a — e, h denoting whether an electron or a 
hole travels on the trajectory, and a subscript identify- 
ing whether the segment touches a superconducting lead 
(a = Si) or only normal leads (o = n). With this conven- 
tion, the blue trajectory in Fig. [1] represents trajectories 
coded in Eq. ([5]). At e = 0, the action phase accumulated 
along such trajectories is given by Sjip) = p{—TT±(l)) + (p, 
where ip is a constant phase that is irrelevant for trans- 
port. The crucial element is that the phase difference 
between the two superconductors cancels out the accu- 
mulated Andreev scattering phase when tt = (f>. Then all 
members of a family interfere constructively with each 
other, because Sj{p) — Sj{p') = 0. This holds simul- 
taneously for all families of trajectories that touch both 
superconductors, with the topology of the trajectories 
sketched in Fig.[T][i]. As there are infinitely many such 
trajectories, the result is macroscopic resonant tunneling 
with a conductance G{(t> — n) (x GoNn, independent of 
r„ ^ 1. Neither macroscopic resonant tunneling, nor 
the associated massive quasi-degeneracy of energy levels 
of Andreev billiards around e = for </> = tt have been 
noticed in earlier investigations of the density of states of 
Andreev billiards 0, 0, H 0, and of transport through 
the Josephson junction set-up of Fig. [T] [ll. 

We sketch our analytical calculation. In the symmetric 
configuration we consider, where each normal lead car- 
ries the same number of channels, connected to the cav- 
ity with the same transparency, the average conductance 
from L to R reads, to leading order in Fn^n uA\ 



{G)/Go = (T- ) + 



LL/ 



(2) 



Here, gives the transmission probability for a quasi- 
particle of type a from the normal lead « to a quasi- 
particle of type /3 into the normal lead j. To evaluate the 
resonant contributions to (T^j^) and (T^£), we follow the 
semiclassical approach of Ref. [l^ (see also Refs. [lij). 
We first write the transmission probabilities as 



(3) 



71,72 



This expression sums over all classical trajectories 7I and 
72 entering the cavity at yo on a cross-section of lead i 
and exiting at ?/ on a cross-section of lead j, while con- 
verting an a quasiparticle into a (3 quasiparticle. The 
phase SS = S-^i — 8^2 gives the difference in action phase 
accumulated along 7I and 72. In presence of tunnel bar- 
riers, the stability amplitude A~f is given by [13, 



A^^B^t,tj W[rk\ 



where l^ (fc) gives the number of times that 7 is reflected 
back into the system from the tunnel barrier k = L, R, 
the transmission and reflection amplitudes at the normal 
leads satisfy ji^p = (1— = ^n,s (for i = L, R, Si, or 
S2), and = {dpyg/dy)^ measures the rate of change 
of the initial momentum py^ as the exit position y of 
7 is changed, for a fixed sequence of transmissions and 
reflections at the tunnel barriers. 

We use Eqs. ^ and (jl]) to evaluate the contributions to 
the total conductance, Eq. ([2|), arising from trajectories 
touching both superconductors such as those sketched 
in Fig. [1] These are the only trajectories that are 0- 
depcndent, they are subdivided into class I trajectories, 
contributing to (T|j£) (blue trajectory on FigU]), and 
class II trajectories, contributing to (T^x) i^'^'^ trajec- 
tory on FigU]). From now on, we focus our discussion on 
class I trajectories. The calculation of class II contribu- 
tions proceeds along the same lines and will be presented 
elsewhere . Class I trajectories are made of the fol- 
lowing sequence 



Ap) _ (e) Jh] [{h) (e) (e) (h) 

7/ — Isl + 7sl + P >^ 7s3 + /s3 + 'si + I si 



(5) 



(4) 



where si and s3 can be interchanged, and p = 0, 1, 2, . . . 
They undergo 2p + 1 Andreev reflections, 2p reflections 
at tunnel barriers [H], and accumulate an action phase 

S-yj = p{-Tr-(j) + etej) + 2ety^,-{TT/2 + (j)/2).{6) 

One should substitute (/> ^ —(p when interchanging seg- 
ments si and s3, but the relative sign between tt and (j) 
does not affect the final result. Here, i^j gives the dura- 
tion of the Andreev loop [the sequence between bracket 
in Eq. ([5])], i^^^ the duration of the segment 7^^. We 
see that at e = and (p = n, the phase difference accu- 
mulated by any two members (with different p and p') 
of a given family vanishes, so that all pairs of trajecto- 
ries within a given family resonate. There is however no 
resonance between members of different families. 

In normal chaotic billiards, the stability of periodic 
orbits decreases exponentially with the number of times 
the orbit is traveled [2^ . The situation is fundamentally 
different in presence of superconductivity, where Andreev 
reflections refocus the dynamics. The stability of a tra- 
jectory is then given by the product of the stabilities 
along the primitive segments (7^1 and 75,3 for class I, 7^1, 
7„ and ^s2 for class II) that the trajectories are made of, 
independent of p . This is true as long as half the du- 
ration of the Andreev loop is shorter than the Ehrenfest 
time te, i.e. the time beyond which an initially narrow 
wavepacket can no longer fit inside a superconducting 
lead [1, [§] . For a quantum dot of linear size Lc and Lya- 
punov exponent A (in absence of superconductivity) , one 
has Te = ln[iVg^/A:FLc]; which determines the rela- 
tive measure of trajectories contributing to macroscopic 
resonant tunneling, together with the average time td 
between two consecutive Andreev reflections. 
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We are now ready to evaluate the dominant contribu- 
tions to conductance close to resonance at e = arising 
from class I trajectories. We start from Eq. ([3]), and, 
following the above considerations, we substitute 



71,72 



(7) 



71,72 



7— prim 

To obtain ([7|), we paired trajectories by class, noting that 
for a given class, all trajectories have the same stability 
but differ only by the number of Andreev reflections at 
the superconductors and normal reflections at the normal 
leads, and by the different action phases they accumulate 
along their Andreev loop. The sum over classes is then 
represented by a sum over primitive trajectories, and the 
exponents a = 1 and c = 1 for class I are determined by 
the number of Andreev and normal reflections in Eq. ([5]). 
Reflection phases at the tunnel barriers do not appear 
because all trajectories are traveled as many times by an 
electron as by a hole. The evaluation of J2 ^7 proceeds 
along the lines of Ref . [l^ , and details will be presented 
elsewhere [l^. The resonant part of the conductance 
from class I and II contributions finally reads 



/rphe \ _ ^^n^n 



(8a) 



X ( 1 - (1 + te/td) exp[-TE/rD] 



/rpee \ " n ^ *n 



1 - 2 r, (1 - r„) cos[7r - 0] + r2 (i - r„)2 ■ 

(8b) 



8N, V2r„A^n + 2r,iv, 

X (1 - (1 + te/td + T|/2Tg) exp[-TE/TD] 

^ 1 + rg (1 - rn)2 

1 - 2 r, (1 - r„)2 cos[^ - 0] + r2 (i - r„)4 ' 



The sum of Eqs. ((8a[) and (j8bp gives the dominant semi- 
classical contribution to the conductance. It exhibits the 
functional dependence of resonant tunneling [compare to 
Eq. ([T])], where the resonance is however always at the 
Fermi level, and is achieved by setting the phase differ- 
ence between the two superconductors at </> = tt. This 
resonance condition is the same for all trajectories. This 
is why the resonance is macroscopic, (x and not of 
order one, as is the case for standard resonant tunneling 
in chaotic systems. In most instances, (T^£)r ^ (^i?L)r- 
Then the resonance height at large te/td, small and 
Eg = 1 is given by G{tt) ~ 7riV„/16. Simultaneously, the 
sharpness of the resonance peak, measured by its width 
at half height, is proportional to En. We also note that 
the effect disappears if the superconductors are poorly 
connected to the normal cavity. Eg — s- 0, as should be. 



The conductance is the sum of the semiclassical con- 
tributions, Eqs. ([8]), and of quantum universal contribu- 
tions. We calculated the latter using Nazarov's circuit 
theory ^ and obtained G„ct(7r)/Go = EniVn/2 [l|. In 
the tunneling regime, this is smaller than the semiclas- 
sical contribution by a factor oc En ^ 1, i.e. semiclas- 
sical contributions enhance the conductance by a factor 

We briefly confirm our predictions numerically. We 
extend the open kicked rotator of Refs. 0, ll^ to take 
into account both transport between two normal leads 
and Andreev reflection at two superconducting terminals. 
We construct a four-terminal scattering matrix from the 
Eloquet operator of the kicked rotator as in Refs. 
and use the method of Ref. |l2l| t o evaluate the exact 
expression for the conductance , 



(G)/Go = rj 



RL 



T. 



he 
RL 



rphe rphe rphe rphe 



'' rphe 



rphe 
^RR, 



rphe 
^LR 



(9) 



In our numerics, we restrict ourselves to perfectly con- 
nected superconductors. Eg = 1. We average our data 
over ensembles of systems with fixed classical parameters 
- such as the width of the leads, the strength of the tun- 
nel barriers and the size and chaoticity of the cavity - but 
different chemical potentials or lead positions. We focus 
our investigations on the semiclassical limit k-pL^ 00. 

In Fig. [21 we show a resonance in the semiclassical 
regime. We obtain very good agreement between the nu- 
merical data (circles) and the analytical prediction (red 
solid line) with te/td — 0.79. Without the semiclassi- 
cal contribution, this agreement would break down close 
to resonance, where universal contributions give a pre- 
diction Gnct(7r) = 2 (green line), too small by an or- 
der of magnitude. The left inset illustrates the increase 
of the peak height and narrowness as the semiclassical 
parameter fcpLc increases, all classical parameters being 
kept constant. The four sets of data in this inset cor- 
respond to a given classical configuration, with the elec- 
tronic wavelength decreasing by factors of four from one 
curve to the next, starting from the bottommost (black) 
curve. The conductance increases at each step because 
the number of conduction channels scales linearly with 
kpLc- In absence of semiclassical contributions, these 
four curves would exhibit the same peak-to-valley ratio, 
but here they obviously do not. This is quantified in the 
right inset to Fig.[2l where we show both the peak height 
and the peak-to-valley ratio corresponding to the same 
configuration as in the main plot, while varying kpLc- 

The connection can be made between the predicted 
and observed enhancement of conductance at ^ = tt 
and resonant tunneling through a macroscopic number 
of quasi-degenerate Andreev levels. In 7r-biased closed 
chaotic Andreev billiards, Bohr-Sommerfeld quantization 
predicts that all periodic orbits touching both supercon- 
ductors contribute to a peak in the density of states at 
the Fermi energy with c>c Ng states. Once electrodes are 
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G{(f>) G(tv) G(7r)/G(0) 



G(0) 




cI>/2tt 

Figure 2: (Color online) Conductance through a chaotic 
Josephson junction vs. the phase difference </> between 
the two superconductors. Circles are numerical results 
obtained from the kicked Josephson rotator with Fn — 0.01, 
fcpLc = 8192, kpLc/Nn = 20, kpLc/Ns = 10 and Lyapunov 
exponent A ~ 1.3. The red curve is the analytical prediction 
obtained by summing the semiclassical resonant contribu- 
tions of Eqs. with the universal prediction. The green 
line gives the universal prediction obtained from circuit 
theory. Left inset; Numerical data for the same classical 
parameters kpLc/Nn ~ 20, kpL^/Ns — 10 and _fL = 10 as 
in the main plot, for kpLc — 128 (black curve), 512 (red 
curve), 2048 (green curve) and 8192 (blue curve). Note the 
change in peak-to-valley ratio. Right inset: peak-to-valley 
ratio G{n) /G{0) (black circles) and peak conductance G(7r) 
(red squares) as a function of kpLc, for the same classical 
configuration as in the main plot. Data are averaged over 
150-1000 sample realizations. 

connected to the billiard, all those e = each level that 
significantly overlaps with the electrodes contributes one 
perfect conductance channel to transport via resonant 
tunneling, which therefore becomes macroscopic. We 
have numerically checked that the observed increase of 
conductance is accompanied by the emergence of a large 
peak around e = in the corresponding Andreev billiard. 
This and other results will be presented elsewhere [l^. 

In summary, we investigated semiclassically the con- 
ductance through quantum chaotic Josephson junctions 
connected to two external normal leads. We found 
an order-of-magnitude enhancement of the conductance 
when the two superconductors have a phase difference of 
TT. We identified the mechanism behind this enhancement 
as resonant tunneling through a macroscopic number of 
quasi-degenerate levels at the Fermi energy. 
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